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Abstract 

An extended version of the pseudo-SU(3) model which includes both spin and 
proton-neutron degrees of freedom is used to study the influence of the pairing 
interaction on K-band mixing, B(E2) values and quadrupole moments. Using the 
asymmetric rotor model as a backdrop, specific consequences of a many-particle 
shell-model based description of these collective properties are demonstrated and 
fundamental limits of the collective model's approach are investigated. Finally, the 
pseudo-SU(3) model, including representation mixing induced by pairing, is used to 
calculate the energies of 140 Ce and the results are compared to experimental data 
and other theories. 



1 Introduction 

Since the discovery of the pseudo-spin symmetry 25 years ago jl], it has attracted the 
attention of numerous physicists and led to many successful applications of the theory 
(see Ref. || for a review). The pseudo-SU(3) model j| takes full advantage of the pseudo- 
spin symmetry and has been used in the description of a wide spectrum of nuclear physics 
phenomena, ranging from collective excitation spectra || and the scissors mode ||, to 
identical superdeformed bands 0, and most recently to /3/3-decay ||. 

In spite of these successful applications, most calculations that have used the pseudo- 
SU(3) model are very schematic in the sense that: A) the Hamiltonian has a rather 
simple structure as it usually includes only a long-range quadrupole-quadrupole term plus 
a residual interaction designed to accommodate certain special properties of low-energy 
collective spectra, and B) the configuration space was severely truncated, usually down 
to only one irreducible representation (irrep) of SU(3). The reason for invoking these 
simplifications was technical, in particular, it was not possible to calculate the matrix 
elements between states of different SU(3) irreps. Very recently a code was released 
that lifts this limitation ||. It allows: A) the inclusion of short-range interactions, like 
pairing correlations, which other shell-model theories indicate to be essential, and B) 
larger configuration spaces which include many SU(3) irreps. 



As a consequence of these developments, in a recently published paper [1C] an extended 



version of the pseudo-SU(3) model, which takes the spin degrees of freedom in a full 
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proton-neutron formalism into account, was reported. Specifically, the pairing interaction 
was expressed in terms of SU(3) tensor operators and the effects of pairing correlations 
on the low-energy collective spectra and moments of inertia were demonstrated. 

It is the purpose of this contribution to extend these studies by investigating the influ- 
ence of pairing on i^-band mixing, B(E2) values, and quadrupole moments. These studies 
will be done by comparing the pseudo-SU(3) model results to those of the asymmetric ro- 
tor model. The latter is a simple version of the geometric collective model which explains 
low-energy nuclear properties in terms of rotations and surface vibrations of a liquid-drop 
type nucleus and neglects all single-particle degrees of freedom. As a phenomenological 
model it has been very successful in the description of a variety of low-energy properties 
of even-even nuclei [[nj. By comparing results of the pseudo-SU(3) and the asymmetric 
rotor models, it is possible to identify fundamental limitations of the collective model 
approach on the one hand, while demonstrating the ability of the pseudo-SU(3) model 
to describe collective nuclear properties correctly on the other. These schematic studies 
provide an understanding of the intrinsic properties of the pseudo-SU(3) model that is the 
basis upon which extended application of the theory for describing and predicting nuclear 
structure phenomena will rest. As a first attempt in this direction, the excitation energies 
of the semi-magic nucleus 140 Ce are calculated and compared to experimental data and 
the corresponding results for other theories. 

This paper is organized as follows: In the next section a brief review of the pseudo- 
SU(3) model and the asymmetric rotor model is given, including a listing of the wave- 
functions, Hamiltonians, and transition operators for each theory. In Section 3.1, the 
importance of K, the intrinsic z-axis angular momentum projection is discussed, and how 
states of different K values are coupled in both the asymmetric rotor model and in the 
pseudo-SU(3) model is explored. In Sections 3.2 the dependence of B(E2) values and 
quadrupole moments on the pairing strength is investigated; specifically, these measures 
as revealed by the pseudo-SU(3) model are compared to the corresponding asymmetric 
rotor model results. In Section 4 the pseudo-SU(3) model is used to describe the ex- 
perimental energy spectrum of 140 Ce and the results are compared to other theories. A 
summary and conclusion is given in Section 5. 



2 Models and Observables 

An outline of the assumptions and definitions of the wave functions and Hamiltonian of 
the pseudo-SU(3) model (Subsection 2.1) and asymmetric rotor model (Subsection 2.2) 
are given in this section. Since the behavior of B(E2) rates and quadrupole moments are 
presented in Section |3], definitions of these observables are given in Subsection 2.3. 

2.1 Pseudo-SU(3) model 

In a recent publication, an extended version of the pseudo-SU(3) model was introduced. 
The formulation that was introduced included general expressions for basis states and for 
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matrix elements of generic operators between them JlO| . Therefore, the present discussion 



is restricted to a brief summary of the essential ingredients of the model. The interested 
reader can find additional details in Refs. 0, |J. 

The extended version of the pseudo-SU(3) model introduced in Ref. |T(J and which 



is employed here, explicitly includes the spin degree of freedom and the Pauli Exclusion 
Principle in a full proton-neutron fermion formalism. Specifically, the pseudo-SU(3) model 
is a Otuv theory; vertical couplings to higher shells are only allowed in an extension of 
the pseudo-SU(3) model called the pseudo-symplectic model ||. In this contribution 
the protons (n) and neutrons (z/) occupy the real-space shells N n = 4 and N u = 5 or, 
equivalently, the pseudo-space shells N n = 3 and N v = 4 which is characteristic of rare- 
earth nuclei. 

The wavefunctions of the pseudo-SU(3) model are classified by Casimir invariants of 
the group chain (see Ref. JED!) 



U(2Sl„) x U(2Q V ) d \u(Q n ) x £/(2) w l x \u(Q u ) x U{2) 



D 



SU(3) {7v) x SU(2)M] x \sU(3) (u) x SU{2) {u) 



D 



S£/(3) w x Sf/(3) M 1 x \SU(2)W x SU(2)M 



D SU(3) x SU(2) S D SO(3) L x SU{2) S D SU(2\ 



where 2Vt a denotes the total number of single-particle levels in the N a shell for protons 
(a = 7i") and neutrons (a = u), respectively. The irreducible representations (irreps) of 
U(2Q a ), labeled by [l m<T ], where m a denotes the number of a-type particles, are the anti- 
symmetrized many- particle proton or neutron wavefunctions which separate into spatial 
(— ► U(Q a ), irrep label [f a ] ) and spin (— > f/(2)( f7 \ irrep label f a ) degrees of freedom. 
All calculations presented in this contribution are restricted to the most symmetric and 
energetically lowest irrep of U(Q a ), which means that couplings to (S ^ 0) modes are 
excluded. This seems reasonable since those spin-flip excitations are significantly higher 
in energy than the low-energy excitations which are of primary interest in this contribu- 
tion. All SU(3) irreps (X a , fi a ) that are contained in a given [f a ] are determined from the 
reduction U(VL a ) D SU(3)^ where possible multiple occurrences of the same (\ a , fx a ) are 



labeled by the integer a a . The reduction 
SU(3) irreps of the total wavefunction, labe 



SU(3)W x S77(3)M -> SU(3) leads to the 
ed by (A, /i), which are determined by taking 



all possible products {(X^, fin) x fi u )} — > p (A, fi) into account, where the running 
index p = 1, . . . , p ma x numbers possible multiplicities. The orbital angular momentum 
values L which are contained in a fixed (A, p)-irrep are determined through the reduction 
SU(3) D SO(3)l and multiple occurrences are numbered by the index k. The total spin 
S is derived from the SU(2) product S n x S v — > S where proton and neutron spins 
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and S v label the irreps of SU(2) n and SU(2) U respectively. Finally the total angular 
momentum J is obtained through the SU(2) coupling L x S —>■ J as symbolized by the 
group reduction [SO(3) L x SU(2) S ] -> SU(2)j. 

The eigenvalues of the Casimir invariants provided by this group reduction scheme, 
augmented with the necessary multiplicity labels, allows for a unique classification of the 
pseudo-SU(3) model wavef unctions, 

\{m w [fir\air{K, /iw) ,m v [f v ]a„(\ v , fi v )}p(X, h)kL {SV, S v } S;JM) 

= XX (^t> A*7r) k-xL % M Lv ; (A M Lu | (A, //) kLM l ) p 

{-} 

x ( ^M<^, S„M & | SM S ) x ( LMl, SM s | JM) 

M Lv S y M Sv ) 

(2) 

with the abbreviation {-} = {M 57r , Ms„, M S} M Ln} M Lv , M L , k^, k v , L w , L u } . In this 
result (LxMl-l-, L 2 Ml 2 \ LMl) denotes SU(2) Clebsch-Gordan coefficients and the SU(3) 



coupling coefficient (. . . ; . . . | . . .) is defined in Ref. [12 



2.1.1 Hamiltonian and Transition Operators 

The pseudo-SU(3) model Hamiltonian, Hpsu(3), used for both the schematic calculations 
in Section ||| and a comparison to experimental data in Section is a many-particle 
extension of the Nilsson-model Hamiltonian, 

Hpsu® = H -^Q a -Q a - G n Hp - G U H" P + aK 2 + bJ 2 + D n ]T l\ + D v £ £ (3) 

where Q a denotes the algebraic quadrupole operator and \ t ne quadrupole-quadrupole 
interaction strength parameter. (Hq denotes the spherical harmonic oscillator which, 
however, is trivial since the configuration space of the pseudo-SU(3) model is restricted 
to a single shell in both proton and neutron spaces. Since only excitation energies are 
of interest, the Hq contribution, which is equal for every configuration, is of no physical 
importance.) Hp (Hp) stands for the pairing interaction for protons (neutrons) which 
is multiplied by the pairing strength parameter G n (G u ). The term K 2 is introduced to 
accommodate the i^-band splitting observed in the low-energy spectra of heavy deformed 
nuclei and J 2 is the square of the total angular momentum. These five terms constitute the 
Hamiltonian for the schematic calculations in Section |3] and have already been discussed 
to some extend in previous publications so the reader is referred to Refs. [[], [H| [TJJ for 



additional details. The parameter set is identical to the one used in Ref. |T(J, namely, 
X = 4.32 keV, a = 202 keV, and b = 9.26 keV, D n = D v = and is compatible with sets 



obtained from best-fit calculations for rare earth nuclei [15 . 
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In the Nilsson Hamiltonian the squared single-particle angular momentum operator, 
I 2 , is of great importance since it effectively flattens the harmonic oscillator potential for 
states of higher angular momentum and in so doing mimics a radial potential shape of 
the Woods-Saxon type. It is therefore necessary to add analogous terms to Hpsu(3) if the 
pseudo-SU(3) model is used for the description and prediction of experimental data. This 
is why the Hamiltonian of Ref. [[HJ has been extended to include the terms D n Yi w l 2 n 
and D v YU* If for protons and neutrons, respectively. The particular choice of parameters 
used for the pseudo-SU(3) model application to 140 Ce is discussed in Section |J 

An expression for the electric quadrupole transition operator and the magnetic dipole 
operator for the pseudo-SU(3) model is given in Ref. @. The result for the real-space 
electric quadrupole transition operator is 



T 2fl (E2) = bl UY.rlmMM) + z«Y;rl($Y 2 »W^ 1 (4) 

\ ien lev / 

where 60 — Aafm is the harmonic oscillator size parameter, e a is the effective nuclear 
charge for protons (cr = n) and neutrons (a = v>), respectively (see Section the spherical 
single-particle coordinates of the valence nucleons are denoted by r a (i),i9f, and (p°, and 
Yin stands for a spherical harmonic PB| . Since all calculations of the pseudo-SU(3) model 
are performed in the pseudo-space, this real-space operator must be expressed in terms 
of pseudo- space quantities. The explicit procedure and the results of this real-space — > 
pseudo-space transformation are given in Ref. 0. 

2.2 Asymmetric Rotor Model 

For the sake of completeness, the most important definitions of the asymmetric rotor 
model are given next. The reader is referred to Refs. [17, ITR [TJ| ^TJ for a more 



detailed account. The asymmetric rotor model Hamiltonian Harm is defined as 

3 T 2 



Harm - ~~T 
k=l Jk 



1 \ I 2 

1 \ , J 3 



where l\ denotes the intrinsic /c-component of the total angular momentum / and /+ (/_) 
the corresponding raising (lowering) operators. The hydrodynamic moment of inertia 
Jk = 4:Bp 2 sin (37 — ^p) depends on the mass parameter B and the quadrupole deforma- 



tion variables 7) (with 0° < 7 < 60°, j3 > 0, see Ref. [F2TJ) all of which are treated as 
adjustable model parameters. The eigenfunctions of Harm, ^ ARM ($i), are determined 
by diagonalization within an orthonormal basis of symmetrized wave functions 
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where the Wigner functions ^[^(i?,) |22|] depend on the three Euler angles (#1, $2, $3) 
and M (K) denotes the laboratory (intrinsic) angular momentum z-axis projection. For 
a fixed angular momentum / the basis dimensionality is given by the number of possible 
K values: K = 0, 2, 4 J if J even and K = 2, 4 7 — 1 if 7 odd. 



Figure [j] 



Finally, recall that the asymmetric rotor model eigenenergies have a characteristic 7 de- 
pendence (see Fig. [I]) which is symmetric around 7 = 30°. (The other two model param- 
eters (3 and B induce only an overall scaling of the eigenenergies.) 

The quadrupole operator Qf? M of the ARM is given in lowest order by 

C = 4* feo(^) cos 7 + (D%(A) + DlUm) ^0 (7) 

In this expression the normalization coefficient A = 3ZRl/4ir depends upon the charge 
Z and the nuclear radius R Q . 



2.3 Physical Observables 

For the sake of completeness, and before continuing with a more detailed consideration of 
the pseudo-SU(3) wave functions, the definitions of some physically important quantities 
will be given. Reduced matrix elements of the SO (3) tensor operator, between states 
with initial (final) angular momentum and projection Jj and Mi (J/ and Mf) are defined 



by 122 



(l}Jf\\Q\\\liJi){JiM h \n\JfM f ) = (jfJfMflQx^iJiMi) 

where 7« and 7/ represent additional quantum numbers that are required to uniquely define 
the initial and final state, and (lfJj\\Q\\\liJi) stands for the reduced matrix element. 
(The use of J for the total angular momentum is customary in microscopic work while 
/ is normally used in macroscopic theories. Here and in what follows, J and I will be 
used interchangeably for the total angular momentum.) The definition of the reduced 
transition probability for electric quadrupole radiation pl| is then given as 



2 7 +1 

B(E2; 7^ - lf J f ) = -L^-^JfWT^W^) 2 (8) 
and the definition of the electric quadrupole moment is 




Qin J) = y ~^~\ n , V 1)(2J ; 3) (7J|OT|| 7 J). (9) 
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3 Pairing and Observables 



The consequences of the pairing interaction within the pseudo-SU(3) model is explored 
in this section by considering its effect on i^-band mixing and on B(E2) values and 
quadrupole moments. The f^-band mixing is considered first because the information 
provided by this effect is easier to realize and simpler to understand than the corresponding 
B(E2) and quadrupole moment results. 

3.1 K-band Mixing 

The quantum number Kj is the projection of the total angular momentum J on the 
system's intrinsic z-axis. Since only S — states are considered in this contribution, 
J = L and Kj = Kl and thus the symbol K can and will be used for both. The K label 
is important in the classification of low- energy collective bands of prolate or near-prolate 
even-even nuclei. Specifically, the ground-state band of these nuclei is associated with 
K = 0, the so-called one-phonon 7 band with K = 2, likewise the two-phonon 7-band 
with K = 4 or K = 0, and so on.[| 

The quantum number K is exactly conserved for prolate nuclei only; couplings between 
bands with different K values are connected with deviations from axial symmetry, that 



is, with 7 7^ 0. From the point of the view of the Geometric Collective Model JIT] , Pq| , the 
.fT-band mixing can be understood within the framework of the very simple asymmetric 
rotor model picture: The Hamiltonian in Eq. |5| is diagonal so long as 7 = 0, since in this 
case J\ = J2 and K is trivially conserved. Increasing 7 (introducing triaxiality) gener- 
ates non- vanishing matrix elements between basis states with AK = ±2. For example, 
the K — eigenstate of the symmetric rotor gains non-vanishing contributions from the 
K = 2 basis function, and so on. This /^-independent phenomenon is illustrated quanti- 
tatively in Fig. |] which depicts the expectation value of K 2 for the even J asymmetric 
rotor model yrast states, (ARM\K 2 \ARM) y , as a function of 7. 



Figure ^ 



The results show that (ARM\K 2 \ARM) y increases with both higher angular momentum 
and larger 7 values: For a fixed value of the 7 deformation, (ARM\K 2 \ARM) y increases 
with increasing angular momentum because the higher the angular momentum, the higher 
the K values (K < I) of basis functions \IMK > that can contribute to the eigenfunc- 
tion. And for a fixed value of the angular momentum, (ARM\K 2 \ARM) y increases with 
increasing 7- value because the non-diagonal term in Harm becomes more and more domi- 
nant. On the far right-hand-side of Fig. |2|, the expectation value of K 2 is indicated for the 
case of complete mixing, that is, if all the diagonalization coefficients have the same ^= 

*Recently the K=A band received some attention through a systematic measurement of 7-vibrational 
anharmonicities in the rare-earth region p3|, B4j 
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magnitude and thus (ARM\K 2 \ARM) y = k J2k A 2 , where d denotes the dimensionality. 
At 7 = 60° these estimates are about 25 percent below the results of the diagonalization, 
indicating that the yrast states favor basis state with higher A values than that of a 
uniform distribution. 

It is now interesting to use these simple results as a backdrop for a study of pairing 
induced A-band mixing in the pseudo-SU(3) model. To accomplish this it is necessary 
to have a shell- model expression for the K 2 operator |T3| , |H] , |26| . For completeness the 
argument that leads to such a result is repeated here. The basic idea follows from a 
frame-independent shell- model expression for the asymmetric rotor model Hamiltonian 
(see Eq. H). To achieve this one can introduce two SO (3) scalars 



A3 = L iQij L j — X] ^ 



2 

i i 



1,3 



XI = Y,L t Q^ k L k = ^X 2 I 2 



where L, and Q?- (A») denote, respectively, projections of the total angular momentum 
and the algebraic quadrupole tensor in either the laboratory or the intrinsic body-fixed 
principal axis system [C|]. (The sums in these expression insure rotational invariance 
and thereby the frame- independent character of the operators.) Adding the total angular 
momentum L 2 = I 2 = J2i If to Ag and A4 means a mapping onto the asymmetric rotor 
model Hamiltonian can be established: 



+ 



1 2 



+ 



7-2 
J 3 



2j7i 2 J/2 IJz 



ah 2 + bX* + cXt 



where the coefficients a, b, and c are real numbers. The geometrical expression for the 
K 2 = I 2 operator emerges trivially from this expression by setting \j1J\ = 1/2 J 2 = 
and 1/2 Jj, = 1. Taking advantage of this observation, it is straightforward to derive a 
shell-model expression for K 2 



T3f 



K 



X1X2L 2 + A3A3 + A4 
2A| + A1A2 



(10) 



where the principle axis quadrupole components can be shown to be given by [|T^] 



Ai = -(-A + /i) , A 2 



I(_ A -2//-3), A 3 = i(2A + ^ + 3). 



The matrix elements of K 2 in the extended version of the pseudo-SU(3) model which is 
used in this contribution can be easily derived using the generic formulas in [14J . 

To investigate the influence of the pairing force strength on A-band splitting in 
the pseudo-SU(3) model, consider a configurations with two protons and two neutrons 
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(27r, 2v) in the real [pseudo] space shells (N n = 4, N v = 5) [(N n — 3, N v — 4)]. This con- 
figuration exhibits a rich structure and its dimensionalities are small enough that basis 
space truncation measures need not be invoked, that is, all possible couplings of proton 
and neutron 577(3) irreps can be taken into account for all angular momentum values. 

Figure |3| 

Figure || shows the expectation value of K 2 for even angular momentum yrast states 
of the (2-7T, 2v) configuration, (27r2z/ | K 2 | 2n2i>) y , as a function of the pairing strength 
parameter G, where for simplicity G n and G v were set equal (G n = G v = G) in the 
calculations. Before entering into a detailed discussion of these results, note that only the 
pairing terms in Eq. [| couple states belonging to different SU(3) irreps (A,//). If G = 
the quadrupole-quadrupole interaction Q a ■ Q a ~ 4C*2 — 3L 2 dominates and forces states 
of the SU(3) irrep with largest second order Casimir invariant G% (Eq. |12|) to lie ener- 
getically lowest. For the (2tt, 2v) configuration this leading irrep is (A,//) = (14,0) which 
contains K = states only and is the reason why the results show {2-k2v\K 2 \2t[2v) y = 
for all value of the angular momentum (Fig. [3])|] 

In a recent paper (see Figures 3-5 in JED]) the intensity distribution of the yrast eigen- 
states of Hp SU ( 3 ) was shown to spread over more and more basis states with increasing 
pairing strength. As a consequence, the (2ir, 2v) yrast states have contributions from 
(A / 0,/j ^ 0) SU(3) irreps which contain K ^ configurations. This admixture shows 
up (Fig. []) as an increase in the (2tt2v\K 2 \2tt2v) y value with increasing pairing strength 
G. (The L = 10, 12, 14h states are unaffected since they are unique.) 

Another mechanism that acts in the same direction is the coupling to higher K states 
within one (A ^ 0, /i / 0) irrep through the X% and X% operators contained in K 2 
(Eq. |10]). This effect increases the expectation value of K 2 for states of higher angular 
momentum as long as there are states with higher K values available. 

For fixed G the (2n2v\K 2 \2^2v) y value increases up to J = 6h and goes down for 
higher angular momenta due to the decreasing availability of higher K basis states. This 
mechanism counteracts the two aforementioned effects and is illustrated in Fig. f| where 
the average K value, K, for a fixed even angular momentum value is given for the (2ir, 2v) 
configuration. 

Figure |] 

Note that the (2ir, 2v) distribution has its maximum at J = 4 and clearly indicates the 
decreasing availability of higher K basis states with increasing J. 

^The selection rule for K in a fixed (A, /Lt)-irrep is given by K = min(A,/i), min(A, /i) -2, min(A,/x) — 
4,...,0or 1 |3 . 
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Figures ||]a) and ||b) 



General results of this type are shown in Figs. [5] where K values for identical particle 
configurations in the N = 3 (Fig. |5|a) and N = 4 (Fig. ^o) shells are given. Only results 
for particle configurations are given since the K values for holes are identical to those 
for particles if the hole and particle numbers are equal. The latter is the reason why 
within a single shell the maximum K value, K max , increases until the shell is half full 
and decreases in a particle-hole symmetric fashion when the shell is more than half full. 
By comparing results for the N = 3 and N = 4 shells one finds that the K distributions 
are shifted towards higher K values with increasing shell number and that as a general 
rule the maximum average K value for each configuration is found at roughly one third of 
the maximum possible angular momentum. This feature is a consequence of the selection 
rule on K values in a fixed 577(3) irrep, which in turn follows from the Pauli Exclusion 



Principle (see Ref. |[28||). One does not expect to find this effect in the phenomenological 
asymmetric rotor model which knows nothing about the single-particle structure. And, 
indeed, there are no rules other than (K < I) which limit the K projection in the 
asymmetric rotor model and thus higher K values can be obtained. 

At the beginning of this section it was noted that a deviation from prolate axial 
symmetry is a necessary condition for the generation of TT-band mixing. This result is 
depicted in Fig. EJ where the 7 dependence of the (ARM\K 2 \ARM) y value is shown. 
It is therefore interesting to consider changes in the 7 deformation that result when the 
pairing strength is increased. The 7 of a particular pseudo-SU(3) model wave function 
can be determined by taking advantage of the mapping between the nuclear deformation 
variables 7) and the S77(3) irrep labels (A, fi) [EM. EUJ: 



[An 1 1 

(11) 

1 -x / C 3 \ 
7 = - cos — 73 » 

where Ci and C3 denote the second and third order Casimir invariants of S77(3) with 
eigenvalues 

C 2 (A, fi) = (A + /x + 3)(A + /x) - A/x, 

(12) 

C 3 (A, //) = (A-/i)(A + 2/i + 3)(2A + /i + 3). 

Since C 2 and C 3 are diagonal in the pseudo-SU(3) model basis functions their expectation 
values are easily calculated. 

Using Eqs. O, the dependence of the 7 deformation on the strength of the pairing 
interaction was determined for the even J yrast states of the (2ir, 2u) configuration. The 
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results are shown in Fig. |6|. 
Figure || 

With no pairing all members of the yrast band which are pure (A, /i) = (14, 0) states with 
the same 7 = 3.2° deformation. But when the strength of the pairing interaction is turned 
on, the system is driven towards triaxiality. As long as the spectrum is dominated by 
the quadrupole- quadrupole interaction (G < O.lMeV) and the rotational character of the 
eigenstates is more or less conserved, the yrast eigenstates are less triaxial the higher their 
angular momentum. Since in this domain the change with increasing angular momentum 
is a relatively small effect, it might be attributable to rotational stretching. For larger G 
values, however, the triaxiality is found to increases more or less uniformly for all of the 
yrast states. For a very large pairing term, where the corresponding 7 values are indicated 
as "Asymptotic Values" on the far right side of Fig. ||, the states show about the same 
7-deformation. 

3.2 Quadrupole Moments and B(E2) Values 

The effect of the pairing interaction on quadrupole moments and B(E2;J — > I + 2) tran- 
sition rates is considered in this subsection. Recall that increasing the pairing strength 
corresponds to increasing the triaxiality of the calculated yrast states of a system. To 
further probe the effect pairing has on collective behavior it is helpful to look at the 7 
dependence of quadrupole moments and B(E2) values. And as above, it is useful to start 
with the complementary asymmetric rotor model results. 

Figure [7] 

Figure § 

Figure [7| depicts the 7 dependence of the spectroscopic quadrupole moments of even J 
yrast states of the asymmetric rotor model. For a range of 7 values induced by increasing 
the pairing strength from weak to strong in the pseudo-SU(3) model , that is, for 7 values 
between about 3° < 7 < 17° (Fig. ||), the negative quadrupole moments are found to 
decrease in magnitude for all values of the angular momentum. Comparing this with the 
corresponding pseudo-SU(3) model results shown in Fig. || shows that the quadrupole 
moments of the even yrast states in the pseudo-SU(3) model also decrease with increasing 
pairing strength. The similarity of these results confirms that an increase in the pairing 
strength goes hand-in-hand with an increase in the 7 deformation. 
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The situation is less clear if the yrast intra-band B(E2) values of the two models are 
compared, see Figures |9] and [10 . 



Figure ^ 
Figure |10 



These figures depict, respectively, the dependence on 7 of the intra-yrast band transitions 
B(E2;J — > / + 2) in the asymmetric rotor and pseudo-SU(3) models. While for 3° < 7 < 
17° the asymmetric rotor model results show a moderate decrease in the B(E2;0 — > 2) 
strength, for all pratical purposes the asymmetric rotor model B(E2;J — > / + 2) values 
are independent of 7. In contrast with this, Fig. [10] shows that the pseudo-SU(3) model 
yields a very pronounced drop in the B(E2;J — > 1 + 2) transition strengths with increasing 
pairing strength. So while diagonal measures suggest that increasing 7 in the asymmetric 
rotor model and the pairing strength in the pseudo-SU(3) model produce similar effects, 
they predict quite different results for intra-yrast band B(E2) transitions. 

4 Application to 5g°Ce82 

The results of the previous section form a backdrop for pseudo-SU(3) model applica- 
tions that describe and predict nuclear properties like excitation energies, B(E2) values, 
quadrupole moments, and ^-factors. A more complete and systematic study of rare- 



earth nuclei will be reported elsewhere this section focuses on the semi-magic H Ce§2 
nucleus which has eight valence protons outside the Z = 50 core. This means the model 
space does not have to be truncated, nevertheless, the 5g°Ce82 system is sufficiently well- 
studied experimentally to allow for a thorough examination of the pseudo-SU(3) model 
description of its properties. 

The 58°Ces2 nucleus has six valence protons in the usual (N n = 4) shell, or equivalently, 
six protons in the pseudo (N^ = 3) space, and no valence neutrons. The neutron part of 
the total wave function is therefore trivial and this in turn implies some simplifications in 
the Hamiltonian of Eq. || Dropping all terms with vanishing matrix elements yields 

Hpsum = ~\Q a -Q a - G V H* P + aK) + bJ 2 +D n J2 1 + ^ (13) 
where the third order Casimir operator, C3, of SU(3) has been added. The physical 



meaning of this operator, see Eqs. [11] and 0, can be easily understood through the 
relation of the pseudo-SU(3) model to the geometric collective model. As mentioned in the 
introduction, this liquid-drop type model describes the low-energy collective excitations 
of even-even nuclei in terms of quadrupole surface vibrations and rotations of the nucleus 
whole. 
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A number of publications show that a phenomenological treatment of the geometric 
collective model provides a very successful means for describing the low-energy proper- 
ties (energies, B(E2) values, and quadrupole moments) of even-even nuclei throughout, 



including phenomena like shape coexistence and shape transition [II], 21]. The basis of 
such an approach is an expansion of the model Hamiltonian in powers of the collective 
((3, 7) variables. A typical ansatz for the nuclear potential in such a scheme is 

V(J3, 7) = x 2 [3 2 + x 3 (3 3 cos 3 7 + x 4 /? 4 + . . . (14) 

where the real numbers X2,x 3 ,X4, ... are usually determined in a least-square fitting 
procedure aimed at an optimal description of the experimental data of the nucleus under 
consideration. 

To understand the physical effect of C 3 in Eq. [13], it is useful to recall the mapping 
between the pseudo-SU(3) labels (A, fi) and the deformation variables (/3, 7) f29| that is 
given in Eq. [11]. The mapping is based on an associated of the invariants of the SU(3) 
and geometric collective models, 

C 2 ~/? 2 , C 3 ~ /? 3 cos 3 7 . (15) 

These relations also provide an explanation for the physical meaning of C 3 in Eq.|i~3|: 
Depending on the positive or negative sign of c in Eq. [T^, or equivalently on the sign of x 3 
in Eq. [TJ], the nuclear system is driven towards more oblate or more prolate deformations. 
So in contrast with the quadrupole-quadrupole interaction, the C 3 term adresses the 7 
degree-of-freedom and helps determine the triaxiality of the system. 

Six parameters (see Table [l]) are required to fix the pseudo-SU(3) model Hamiltonian 
(Eq. ^) for 5g°Ce82- These values were determined in a best-fit calculation based on 
known experimental energies and B(E2) values of the system. 



Table | 



There are a number of theoretical and experimental studies which focused on the 140 Ce 
nucleus |[31| . The current study focuses only on its excitation spectrum; the reader is 



referred to flT5| for other measures as well as a systematic comparison of the properties of 



several different rare earth nuclei. 



Figure 11 



Pseudo-SU(3) model results for 140 Ce energies up to about 3 MeV are shown in Fig. O. 



A comparision with experimental excitation energies that was taken from [32| is given to 
the immediate right of the pseudo-SU(3) results followed by the spectrum obtained by 
Wildenthal |33| and the results of a model study by Waroquier and Heyde |34 . 
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The Wildenthal results were obtained in a conventional shell model calculation that 
used a modified surface 5 interaction in a truncated configuration space. The single- 
particle energies and the parameters of the modified surface 5 interaction were adjusted 
to describe the excitation energies of nuclei in the mass range from A = 136 to A = 145 



33|| . In addition to a surface S interaction, Waroquier and Heyde |34j used a Gaussian 
residual nucleon-nucleon interaction. The single-particle energies were determined by 
solving the inverse-gap equations. 

While the agreement between the three theories is of approximately the same quality 
for energies below about 2.5 MeV, the pseudo-SU(3) model does not account for a low- 
lying 5 + state which experimental evidence suggests lies at 2.349 MeV. Otherwise, both 
the pseudo-SU(3) model and the calculations by Wildenthal account for virtually all levels 
up to about 3 MeV while the results by Waroquier and Heyde seem to miss a few of the 
experimentally established states. In addition, note that the two former models predict a 
few levels, like two 4 + states at about 2.8 MeV, which have not been found experimentally. 
The pseudo-SU(3) spectrum also predicts a + state at about 2.25 MeV that cannot be 
found in the experimental results. 



5 Summary, Conclusion and Outlook 

The effect of pairing on .fT-band mixing, B(E2) values, and quadrupole moments has 
been studied within the framework of a complete (untruncated) pseudo-SU(3) model 
theory. The asymmetric rotor model, which is a collective model theory that describes 
nuclear properties by means of rotations only, was introduced as a backdrop for helping 
to identify collective properties of the pseudo-SU(3) model and, in a complementary way, 
for discovering limitations of the collective model approach: 

• From general considerations it is clear that deviations from axial symmetry are a 
necessary condition for fT-band mixing. This was seen explicitly for the asymmetric 
rotor model by noting that for 7^0 the non-diagonal part of the Hamiltonian is 
non-zero, and for the pseudo-SU(3) model since an increase in the pairing strength 
not only induces K-band mixing but also generates an increase in 7 deformation. 

• The K label of the pseudo-SU(3) model is limited by a SU(3) to SO (3) selection 
rule. This restriction is inherent to the Pauli Exclusion principle and not part of 
the asymmetric rotor model where, in principle, there is no maximum K value. 

• The quadrupole moments of yrast states of the asymmetric rotor model decrease in 
magnitude with increasing 7 deformation (0 < 7 < 30°) and in the pseudo-SU(3) 
model with increasing pairing strength. Since an increase in the pairing strength 
was found to be correlated with an increase in 7 deformation from about 7 = 3° to 
7 = 17°, the physics behind the two phenomena appears to be very similar. 

• The yrast intra-band B(E2;J — > I + 2) strengths decrease rather strongly with 
increasing pairing strength while the asymmetric rotor model values show only a 
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moderate decrease with increasing 7 deformation. 

Finally, the pseudo-SU(3) model was used for a description and prediction of experimental 
energies of 140 Ce and the results were compared to other theories. The Hamiltonian 
parameters were determined in a best-fit calculation that used as input the experimental 
energies and the B(E2)-values of the low-lying states. The pseudo-SU(3) model was found 
to describe the known experimental energies satisfactorily, in a addition, a few new level 
were predicted. A more complete and systematic comparison of the predictions of the 
pseudo-SU(3) model with experimental data and other theories will be reported in a 
forthcoming contribution fT5|| . 

Regarding future developments it should be obvious that the presented formalism can 
be easily applied to (S 7^ 0) states and hence to a description of not only even-even but 
also odd-even, even-odd, and odd-odd nuclei. To put this study into perspective, it is 
also important to note that the results learned from pseudo-SU(3) model studies apply 
as well to the pseudo- symplectic model which is a natural extension that takes couplings 
to higher shell fully into account. 
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JL CAiL Oj1J.1v; Uv/1 


140 Ce 


X [keV] 


3.97 


D„ keV 


-101.05 


G w [keV] 


157.5 


a [keV] 


130.6 


b [keV] 


0.085 


c [keV] 


5.078 



Table 1: List of parameters used in the pseudo-SU(3) calculations of 140 Ce (see text), 
values 
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Figure captions 



Figure 1: The 7 dependence of the lowest asymmetric rotor model eigen- energies in 
units of 4-§^2 . The angular momentum values are indicated on the far right and along the 
central vertical line. 



Figure 2: The 7 dependence of the expectation value (ARM\K 2 \ARM) y of the K 2 op- 
erator is shown for even angular momentum yrast states of the asymmetric rotor model. 
The angular momentum of the various curves are indicated on the far right. The values 
for (ARM\K 2 \ARM) y that are obtained for complete mixing, that is, when all the diag- 
onalization coefficients are set equal to one another, are indicated by horizontal lines at 
the far right labeled above with the corresponding value of the angular momentum. 



Figure 3: The expectation value of the K 2 operator in even angular momentum yrast 
states is shown as a function of the pairing strength parameter G, where for simplicity 
the latter were all set equal, G n = G v = G. The angular momentum value of each curve 
is indicated on the far right. The horizontal lines on the far right indicate asymptotic 
values for the matrix elements (ARM\K 2 \ARM) y that are obtained in the limit of a very 
large pairing strength. 
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Figure 4: This figure illustrates an average property of the pseudo-SU(3) model basis 
states, namely, the average K value for all basis states of a fixed even angular momentum 
J. The K distribution shown is for the (m n = 2, m„ = 2) configuration. 

Figure 5: This figure is similar to Figure 4. The K values were calculated for all even 
angular momentum states J and different particle numbers m in the N = 3 shell [ 5a), 
left] and the N = 4 shell [ 5b), right]. (Particle distributions corresponding to more than 
half- filled shells can be obtained by invoking particle-hole symmetry.) 

Figure 6: The expectation value of the 7 deformation, (7), is shown for even angular 
momentum yrast states of the (m n = 2, m u = 2) configuration. The abscissa is the 
pairing strength parameter G where for simplicity the proton and neutron strengths were 
set equal, G w = G v = G. On the far right the (7) values for a very large value of the 
pairing strength (" Asymptotic Values") are indicated as bars, each of them labeled by 
its angular momentum value. 

Figure 7: Quadrupole moments of yrast states of the asymmetric rotor model for even 
values of the angular momentum are shown as a function of the 7 deformation where the 
numbers on the far left and right denote the angular momentum in units of h. 

Figure 8: Quadrupole moments of yrast states of the (m n = 2, m u = 2) configuration for 
even values of the angular momentum are shown as a function of the pairing strength G. 
The numbers next to the curves on the right denote the value of the angular momentum 
in units of K. For comparison, values for the prolate rotor are indicated on the far left. 
Note that these are almost identical to the pseudo-SU(3) model values. On the far right 
quadrupole moments for the case of very large G are indicated by dashed horizontal lines 
(" Asymptotic Values "). 

Figure 9: The 7 dependence of the intra-band transition probabilities, B(E2;J — > 7 + 2), 
for the asymmetric rotor model yrast states are shown, with the initial and final values 
of the angular momentum indicated on the far right. 

Figure 10: The intra-band reduced transition probabilities, B(E2;J — > I + 2), for eigen- 
functions of the pseudo-SU(3) model are shown as a function of increasing pairing strength. 
On the far right the horizontal lines labeled "Asymptotic Values" are for a very large value 
of the pairing strength. The dotted lines from the smaller G results are included to help 
guide the eye. Similarly, on the far left the horizontal lines labeled "Rotor Values" are 
the B(E2;J — > I + 2) strengths of the prolate rotor. 
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Figure 11: From left to right the figure depicts the excitation spectrum of 140 Ce as 
calculated within the pseudo-SU(3) model , the experimental values (left center), the 
results of a calculation by Wildenthal (right center), and, on the very right, the energies 
as obtained in the model used by Waroquier and Heyde (see text). Note that the thin 
dashed lines between the levels are supplied to guide the eye between corresponding levels 
and that the spacings between the ground and the first excited states are not to scale.) 
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